Topological characterization of quantum phase transitions in a S=l/2 spin model 
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We have introduced a novel Majorana representation of S=l/2 spins using the Jordan- Wigner 
transformation and have shown that a generalized spin model of Kitaev defined on a brick-wall 
lattice is equivalent to a model of non-interacting Majorana fermions with Z2 gauge fields without 
redundant degrees of freedom. The quantum phase transitions of the system at zero temperature 
are found to be of topological type and can be characterized by nonlocal string order parameters. 
In appropriate dual representations, these string order parameters become local order parameters 
and the basic concept of Landau theory of continuous phase transition can be applied. 



The Landau theory of second order phase transitions 
has fertilized modern statistical and condensed matter 
physics. Essential is to use local order parameters to 
describe the continuous phase transition between a dis- 
ordered and an ordered phase associated with symmetry 
breaking, lj However, a quantum phase transition driven 
entirely by quantum fluctuations at zero temperature can 
occur between two disordered phases without any sym- 
metry breaking. 00 A typical example is the topological 
phase transition between two neighboring quantum Hall 
plateaus in the fractional quantum Hall effect. Q As no 
conventional Landau-type order parameters can be used, 
a comprehensive characterization of this kind of quantum 
phase transitions has become one of the most challenging 
issues in condensed matter theory. 

In this work, we present a theoretical analysis of quan- 
tum phase transitions in the following S — 1/2 spin 
model first introduced by Kitaev^ 
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where j and I denote the column and row indices of the 
lattice. This model is defined on a brick wall lattice (Fig. 
IBOj which is an alternative representation of a honey- 
comb lattice (Fig. [TJi). Here we emphasize the use of 
the brick-wall lattice because, as will be shown later, it 
presents an intuitive configuration for labeling sites from 
which both the Jordan- Wigner and duality transforma- 
tions of spins can be more naturally introduced. The one- 
and two- row limits of the brick- wall lattices with periodic 
boundary conditions are simply a single spin chain (Fig. 
ITJi) and a two-leg ladder (Fig. ^p) , respectively. In these 
limits, the analysis of the model is greatly simplified, 
while the results obtained are highly nontrivial and can 
be extended with proper modification to other brick wall 
lattices. In the work of KitaevQ, the two-dimensional 
phase diagram of the ground state was studied. Here we 
will concentrate more on the analysis of quantum phase 
transitions in an arbitrary brick-wall lattice. Our stud- 
ies reveal many generic features of topological quantum 
phase transitions in this system. It also sheds light on 
the understanding of more general systems, where ana- 
lytic solutions are not available. 
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Figure 1: A brick- wall lattice (c) with its equivalent honey- 
comb lattice (d). (a) and (b) are the one- and two-row limits 
of the brick-wall lattices, respectively. 



To solve this model, let us first introduce the following 
Jordan- Wigner transformation to represent spin opera- 
tors along each row with spinless fermion operators: 



+ = Oflt fc<! a lk a ik+E, <J <*«<»«) 
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For each pair of fermion operators (a, a'), we can fur- 
ther define two Majorana fermion operators (c, d): Cji = 



- a,ji) and dji 
aj ; + dji and d 



, + Ojj when j + I 
i{a\l — CLji) when j + I 



even, and 
odd. 



With these definitions, it is straightforward to show that 
at each given row only c-type of Majorana fermions are 
present and the Hamiltonian can be expressed as 

H = -i (JtfjjCj+ij - J 2 Cj-i t iCj : i 

j-\-l—even 

■ -I ;l ),}C,JC ,j . ; ). (3) 

where Dji = idjjdj.1+1 is defined on each vertical bond. 
Since there are no direct connections between any two 
vertical bonds, all Dji are good quantum numbers. Each 
Dji acts like a local static Z2 gauge field. They commute 
with each other and with the Hamiltonian. 

Eq. © is simply a model of free Majorana fermions 
with local Z2 gauge fields. In fact, a similar Majorana 
fermion representation of the Hamiltonian was first used 
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by Kitaev [f| to solve rigorously the ground state in the 
two-dimensional limit. However, in the work of Kitaev, 
a S — 1/2 spin operator is represented by four Majo- 
rana operators, which introduces two redundant degrees 
of freedom at each site and many important properties of 
the system, especially those involving excitation states, 
are blurred by those unphysical states. To remove these 
redundant degrees of freedom, a complicated projection 
on the wave function has to be imposed. 

In Eq.JSJ), Dji for a given I can be all changed into 
—Dji by a unitary transformation. For example, the 
unitary operator U 2 i = Yij c 2j,2i can effectively convert 
all D 2 j : 2i on the (2Z)-th ladder to —D 2 j^i, and U21-1 = 
Yij c 2j-i,2l— 1 converts all £>2j— 1.2/-1 on the (21 — l)-th 
ladder to — D2j-i,2l-t- Therefore, all the eigenstates of 
the Hamiltonian are at least 2 A/ -fold degenerate for a M- 
leg ladder (M = for a single chain). In the discussion 
below, without loss of generality, we will assume all three 
coupling constants J a (a — 1, 2, 3) to be positive. 

For a single spin chain, the Hamiltonian can be read- 
ily diagonalized. The eigen-spectra contain two bands 
of quasiparticle excitations. In the ground state, one 
of the bands is fully occupied while the other is empty. 
When J 1 ^ J2, there is an energy gap in the low- lying 
excitations. By varying J1/J2, we find that the second 
derivative of the ground state energy density Eo diverges 
logarithmically at J\ = J 2 where the gap vanishes (Fig. 

Furthermore, by applying the renormalization group 
analysis to this model, it can be shown that J\/ J 2 will 
flow to infinity if J\ > J 2 or to zero if J\ < J 2 , and 
J1/J2 = 1 is a quantum critical point (Fig. |2t>). 



to rewrite the Hamiltonian £[J as 
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This is nothing but a one-dimensional Ising model with 
a transverse field defined in the dual lattice. When J\ > 
J 2 , it is known that a long range order exists in the dual 
spin correlation function of t 2 j 
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Thus t 2 j can be regarded as an order parameter charac- 
terizing the phase transition from 3\ > J 2 to J% < J 2 in 
the dual space. In the original spin representation, TqT 2 a 
is a string product of <j^: 
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Eq. JBJ then indicates that H has a hidden topological 
order in the J\ > J 2 phase with A x = lim :) _ >00 (A a ;(j)}. 

When Ji < J 2 jthe dual spins t 2 j become disordered 
and A x vanishes 7]. However, by swapping the J 2 - with 
the Ji-term in Eq. (JIJ and applying a similar duality 
argument, it can be shown that a string order of a% 

2j+l 2j+l 

A y = lim (J] &») = (-iy lim< Y[c k ), (8) 
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Figure 2: The second derivative of the ground state energy 
density Eo (a) and the nonlocal string order parameters (b) 
for a single chain. The arrows in (b) indicate the RG flow. 



The transition between the two gapped fcrmionic 
states at J\/J 2 = 1 does not involve any change of 
symmetry. Instead it manifests as a change of topolog- 
ical order. To see this, let us introduce a spin duality 
transformation[£j 

2JV 

^=^-1^ n ) 1R ( 4 ) 

k=j 



is finite in the J\ < J 2 phase and zero otherwise. 

The above discussion indicates that in the single-chain 
limit, the model is in disordered phases, but it contains 
two hidden string order parameters. The quantum phase 
transition at J\/J 2 = 1 corresponds to a continuous 
change of or A y from zero to a finite value from one 
side of the critical point to another. In the dual space, 
however, these nonlocal string order parameters become 
local . This suggests that Landau-type concepts of con- 
tinuous phase transition can still be applied to this simple 
model, but in the dual space. 

For a two-leg spin ladder, the ground state is in a 7r-flux 
phase. Q After diagonalizing the Hamiltonian, wc find 
that the cigcn-spcctra contain four branches of fcrmionic 
quasiparticle excitation bands. At zero temperature, two 
bands are fully occupied and the other two are completely 
empty. In general, there is also an energy gap between 
the ground and excited states. But this gap vanishes 
along both J_ = J3 and J_ = — J3 lines (J± = J\ ± 
J 2 ). As shown in Fig. iJSJd), the second derivative of the 
ground state energy density diverges logarithmically on 
these two lines. 

To clarify the topological nature of these three gapped 
phases as shown in Fig JS^), let us first relabel all the 
sites along a special path as shown in Fig. to express 
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Figure 3: (a)The phase diagram with RG flows for the two-leg 
ladder, (b) d 2 Eo/dJ- (J- = Ji — J 2) along the dotted line in 
(a), (c) and (d) are two deformed single chain representations 
of the two- leg ladder. J3 is set to 1. 




(b) 



Figure 4: (a) The six-row brick wall lattice and (b) its equiv- 
alent three-row lattice. 



the original spin Hamiltonian as an effective single chain 
model with the third nearest neighbor couplings: 



H = T / (J 1 al 
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23+3) 



(9) 



By applying the similar duality transformation @ to this 
Hamiltonian, we find that this model is equivalent to an 
anisotropic XY spin chain with a transverse field in the 
dual space: 



2j) 



(10) 



where W 2 j = T 2j-3 T 2j-i T 2j+i ^ s ^ ne plaquette operator 
defined in the dual space. It is a good quantum number. 

In the ground state, W 2 j = — 1, corresponding to the 
7r-flux phase. It is known that the system is in a long- 
range ordered state of riy when J_ > J3 and in a disor- 
dered state otherwise. [3 Thus, (t^ ) is an effective order 
parameter characterizing the quantum phase transition 
across the critical line J_ = J3. Back to the original spin 
model, there exists a string order parameter A x defined 
by Eq. but in the deformed lattice (Fig. in the 
J_ > J3 phase. Using the results given in Ref. 0, we 
find that 



A, 



2(1 + J+/J_) 



1/4 



(11) 



In the other two phases (J_ < J3), A x = 0. 

The above duality argument can be generalized to clar- 
ify the quantum phase transition across the line J_ = 
— J3 by swapping 3\- with J2 -term in the original Hamil- 
tonian. In the ground state, it can be also shown that the 
nonlocal string order parameter A y defined by Eq.JSJ, 
but now in a deformed lattice where crj is labelled ac- 
cording to Fig. is finite in the phase J_ < — J3, 
and vanishes in the other two phases. Thus the phase 
diagram of the two-leg ladder system can be classified by 



the two nonlocal string order parameters: and A y . 
They are finite in the J_ > J 3 and J_ < — J 3 phases, 
respectively. The quantum phase transition at each crit- 
ical line J_ = ± J3 corresponds to a continuous change of 
one of the string order parameters from zero to a finite 
value. 

For a 2M-row brick wall lattice with M > JL the 
ground state of the model is in a zero-flux phase |8| and 
all Dji for a given / have the same value. The ground 
state is thus 2 2M -fold degenerate, and this massive degen- 
eracy leads to an extensive zero temperature entropy. In 
what follows, we will focus on the case Dji = (— ) . The 
resulting conclusions can be easily generalized to other 
degenerate states by unitary transformations. 

To elucidate the topological nature of the quantum 
criticality in this system, let us first follow what we have 
done for the two-leg ladder to map a 2M-row lattice to 
a M-row lattice by chaining all neighboring sites along 
the stair paths of J\ — J3 links with periodic boundary 
conditions. An example of such mappings is displayed in 
Fig. 0| for a six-row brick wall lattice. In the new lattice 
system, the Hamiltonian 10 can be expressed as 

2N M 

H = — i ( J\C2j-\,gC2j,OL — J2C2j,gC2j+3,a + l 

3 = 1 a=\ 

+ Jzi-Y C2j,aC2j+l,a)- (12) 

By transforming the second index of Cj j0l into its momen- 
tum form Cj, q , Ea. ljT2l then becomes H = J2 q H q , where 



H„ 



■22[J 1 C2j-l- q C2j, q ~ J2e lq C 2j .- q C2j 



+3,9 



+ h{-yc2 .- q c 2j 



+ 1,9.1 > 



(13) 



and q = lirmjM and m = • • • M — 1. This Hamiltonian 
is now block diagonalized according to the value of \q\. 
When q = or tt, — q is equal or equivalent to q, and 
Cj. q is still a Majorana fermion. In other cases, Cj^ q and 
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Figure 5: The phase diagram with quantum critical lines of 
the 2M-row brick-wall lattice. J3 = 1 



Cj t q are conjugate pairs of standard fermion operators, 
i.e. Cj^ q = Cj q . 

The quasiparticle eigen-spectra can be obtained by di- 
agonalizing the Hamiltonian at each \q\ sector (\q\ < tt ). 
Fig. shows the phase diagram of the multi-row lattices 
with M = 2, 3, 4 and 00. For a 2M-row lattice, there are 
M + 2 (M + 1) critical lines if M = even (M = odd). 
These critical lines correspond to the vanishing gap lines 
of quasiparticle excitations, where the second derivatives 
of the ground state energy density still diverge logarith- 
mically. 

When q = 0, H q has exactly the same form as the 
Hamiltonian of a two-leg ladder with a 7r-fiux, repre- 
sented in the deformed lattice in Fig. (J3J;). As discussed 
before, the critical excitation gives rise to two quantum 
phase transitions at the lines J_ = ± J3. The correspond- 
ing topological order parameters, A^o an d Ay 0, can be 
defined from Eas. l|7l8fl by replacing Cj with Cj, q= Q. A X] q 
is finite in the phase J_ > J3 , while A y .o is finite in the 
phase J_ < — J 3 . 



When q = n (M = even only), H q is also equivalent 
to the Hamiltonian of a two-leg ladder, but with a minus 
J2-term. The ground state of H q is effectively in zero- 
flux. Two critical lines are now located at J + = J3 and 
at J_ = with J + > J3. Correspondingly, one can also 
define two topological string order parameters, A X)7r and 



Aj, i7r , from Eqs. J7JSJ) by replacing Cj with c^ q 

(Ay iff ) is finite in the J + > J 3 and J\ > J2 {J\ < J2) 

phase. 

When q =/= or n, Cj^± q are no longer Majorana 
fcrmions. By considering the limiting cases of J\ 3> J 2 
and J\ <C J%, it can be shown that two string order pa- 
rameters defined from the combinations of c^ q and Cj,_ g 
exists. However, as H q is coupled with H- q , we are still 
unable to write down accurately the expressions of the 
string order parameters. 

In the two dimensional limit (M — 00), all the criti- 
cal lines inside the regime enclosed by the three curves, 
J_ = ± J3 and J + = J3 , merge together and form a con- 
tinuum gapless phase. Therefore, the whole phase space 
contains three gapped phases and one gapless phase, in 
agreement with the phase diagram obtained by Kitaev. [j| 
The gapless phase may have a complicated non-Abelian 
topological structure. [3 However, the phase transition 
from any gapped phase to the gapless phase can be char- 
acterized by a string order parameter, since the phase 
boundaries of the gapless phase are fully determined 
by the phase transition lines of the decoupled fcrmionic 
chains with q — and it. 

In summary, using the Jordan- Wigner transformation, 
we have shown that the Hamiltonian defined by Eq. Q i s 
mapped onto a model of free Majorana fermions with lo- 
cal Z2 gauge fields without any redundant degrees of free- 
dom. By solving this model exactly at zero temperature, 
we find that the system undergoes a number of continu- 
ous phase transitions. These quantum phase transitions 
are not induced by spontaneous symmetry breaking and 
there are no conventional Landau-type local order param- 
eters. However, each transition corresponds to a continu- 
ous change of a topological string order parameter from a 
finite value to zero. This reveals that the quantum phase 
transitions in this system can not only be indexed by a 
topological quantum number, but also be characterized 
by a topological order parameter. In the dual space, the 
string order parameter actually corresponds to a local or- 
der parameter and the basic concepts of Landau theory 
of continuous phase transition are still applicable. 

The above conclusion is drawn based on the analysis of 
a specific model l(T)l. however, we believe it might be valid 
generically. Similar topological string order parameters 
have been found, for example, in the Haldane gapped 
spin chains. [jiJ, UJ$ Further exploration of this problem 
is desired, which may lead to a unified description of 
convention and topological quantum phase transitions. 
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